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It is well known that the presence of a weak cross flow in an otherwise two-dimensional 
shear flow results in a spanwise variation in the mean streamwise velocity profile that 
can lead to an amplification of certain three-dimensional disturbances through a kind of 
resonant-interaction mechanism (Goldstein & Wundrow 1994). The spatial evolution of an 
initially linear, finite-growth-rate, instability wave in such a spanwise-varying shear flow is 
considered. The base flow, which is governed by the three-dimensional parabolized Navier— 
Stokes equations, is initiated by imposing a spanwise-periodic cross-flow velocity on an 
otherwise two-dimensional shear flow at some fixed streamwise location. The resulting 
mean-flow distortion initially grows with increasing streamwise distance, reaches a maxi- 
mum and eventually decays through the action of viscosity. This decay, which coincides 
with the viscous spread of the shear layer, means that the local growth rate of the insta- 
bility wave will eventually decrease as the wave propagates downstream. Nonlinear effects 
can then become important within a thin spanwise-modulated critical layer once the local 
instability-wave amplitude and growth rate become sufficiently large and small, respectively. 
The amplitude equation that describes this stage of evolution is shown to be a generaliza- 
tion of the one obtained by Goldstein & Choi (1989) who considered the related problem 
of the interaction of two oblique modes in a two-dimensional shear layer. 
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1. Formulation 


To fix ideas, we consider an incompressible shear flow formed at the interface between 
two parallel streams of differing velocity or alternatively between a single parallel stream 
and a flat plate. The Cartesian coordinate system (x, ?/, z) is attached to the interface with 
x in the direction of the external flow, y normal to the interface, and z in the spanwise 
direction. All lengths are non-dimensionalized by 6+ where 6* characterizes the local shear- 
layer thickness at x = 0. The time t, velocity u — iu + jv -f kw, and pressure variation 
p from the external value P * are non-dimensionalized by U * and p*l 7 * , respectively, 

where £/* characterizes the velocity of the external flow and p* is the density. With this 
non-dimensionalization, the Navier— Stokes equations become 

V-u = 0, (1.1) 

u t + U'Vu + Vp = R~ l V 2 u, (1.2) 

where V = id/dx + jd/dy + kd/dz is the gradient operator, 

R = 8*U*!v*> 1 (1.3) 

is the local Reynolds number, v+ is the kinematic viscosity and an independent variable 
used as a subscript denotes differentiation with respect to that variable. 

The solutions to (1.1) and (1.2) that are of interest here can be represented as the sum 
of a steady base flow plus a time-dependent perturbation, 


u = U(x) + eu(x,t), 

(1.4) 

p = P(x) + cp(x,t), 

(1.5) 
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where e characterizes the local amplitude of the perturbation at x = 0 . Substituting ( 1 . 4 ) 
and (1.5) into ( 1 . 1 ) and (1.2) gives 

VU = 0 , ( 1 . 6 ) 

U-VU + VP = JT 1 V 2 */, ( 1 . 7 ) 

for the base flow and 

V*u = 0, (1.8) 

Ut + U'Vu + tt*V(Z7 + eu ) + Vp = P - 1 V 2 ?2, ( 1 . 9 ) 

for the perturbation. 

The steady spanwise-periodic base flow { U , P} evolves over the long streamwise scale, 

*2 = x/R, (1.10) 

and has an £(<$*) wavelength in the spanwise direction. This implies that the base-flow 
solution expands like 

U = iU 0 (x 2 ,y,z) + R~ l V 0 (x 2 ,y,z) + •••, (1.11) 

P = R~ 2 P 0 (x 2 ,y , *)+••♦, (1.12) 

where V denotes the base-flow velocity in the transverse (or y-z) plane. Substituting ( 1 . 11 ) 
and (1.12) into (1.6) and (1.7) shows that the leading-order base-flow solution is determined 
by the parabolized Navier-Stokes equations (Rudman & Rubin 1968), 

^o ® 2 + V r -V'o = 0, (1.13) 

U 0 (iU 0 + V 0 ) X2 + Vo'V T (iUo + V 0 ) + V T P 0 = V 2 (iP 0 + V 0 ), (1.14) 

where V r = jd/dy -f kd/dz is the gradient operator in the transverse plane. 
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It is assumed that the initial amplitude of the perturbation is small enough so that 
cu < U 0 over the streamwise region of interest. Substituting (1.11) into (1.8) and (1.9) 
then yields 

V-u = 0, (1.15) 

M + t(V r ff 0 -tt) + Vp = 0(R~ l ) (1.16) 

where D = d/di + Uod/dx is the leading-order convective derivative relative to the base 
flow. These equations are just the familiar equations for the linear perturbations about 
a uni-directional transversely sheared base flow (Goldstein 1976; Henningson 1987). It is 
well known that the velocity fluctuations can be eliminated between (1.15) and (1.16) (see 
Goldstein 1976, pp. 6-10 for a detailed derivation) to obtain the following equation for the 
pressure fluctuation 

DV 2 p - 2 V T Uo-V T p x = 0(R~'). (1.17) 

Attention will be restricted to perturbations that are spatially growing and periodic in 
time with, at least initially, a single angular frequency, say T 1 *. The relevant solutions to 
(1.15)— (1.17) then form a spanwise periodic instability wave that propagates in the stream- 
wise direction. The local amplitude of the instability wave increases as the wave propagates 
downstream, but its local growth rate will ultimately decrease owing to the combined effects 
of the viscous spread of the basic shear layer and the viscous decay of the mean streamwise 
vorticity. Nonlinear effects can then become important first within a thin critical layer 
located at the transverse position where the phase speed of the instability wave equals the 
base-flow velocity Uq (once the instability- wave amplitude and growth rate become suffi- 
ciently large and small respectively). In this stage of development, the unsteady flow outside 
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the critical layer remains essentially linear but the instability-wave amplitude is completely 
determined by the nonlinear motion inside the critical layer. 

With this in mind, the origin of the x axis is chosen so that the deviation, 

aSi = S - S 0 < 0, (1.18) 

of the local Strouhal number (or non-dimensional angular frequency) S = 6*F*/U+ from its 
neutral (or zero-growth) value £0 is 0(a) where a <C 1. The precise relationship between e 
and a will be specified below when the flow in the critical layer is analyzed. The relevant 
solutions to (1.16) and (1.17) are then of the form 

u = Re + . . . , (1.19) 

V = ^ (Ape 1 *) + . . . , (1.20) 

where ^4(a?i) is an amplitude function that accounts for the slow growth of the instability 

wave, 

xi = ax (1.21) 

is the streamwise scale over which the wave growth occurs, 

X = a 0 x- St, (1.22) 

is a normalized streamwise coordinate in a reference moving with the wave, and op is the 
neutral wavenumber. The ellipses in (1.19) and (1.20) indicate harmonics of the fundamental 
instability wave that are generated by the critical-layer nonlinearity. Since these harmonics 
do not interact outside the critical layer (to the order of accuracy considered here), their 
outer solutions can be determined a posteriori. 
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Substituting (1.20) into (1.17) shows that, outside the critical layer, the function p of 


X \ , y and 2 is determined to the required order of accuracy by 


' VrP ' 
.(U 0 ~c)\ 


2 ~ 
a p 


(U 0 - cf 


= 0 


(1.23) 


where 


a = Qf 0 — crlA 1 /A, 


(1.24) 


and 


c s S/a, (1.25) 

are the generalized wavenumber and phase speed, respectively, and a prime denotes differ- 
entiation with respect to the argument. It follows from (1.15), (1.16) and (1.19) that the 
velocity fluctuations are determined in terms of p by 


i'(iau) -f = 0, 


(1.26) 


and 


ia(£/o — c)u + i(V T Uo'U + lap) + V T p = 0. 


The solution to (1.23) that satisfies 


p y — 0 at y = 0 ; boundary layer 

p -* 0 as ?/ — » — 00 ; free- shear layer 


and 


(1.27) 


(1.28) 


p — > 0 as y — ► 00 , 


(1.29) 


is analyzed in the following section. 


6 



2. Unsteady flow outside the critical layer 


Outside the critical layer, the shape functions {u,p} expand like 


it = u 0 (y,z) + crui (2.1) 

P = Po(y,*) + crpi(x 1 ,y,z) + •••, (2.2) 


as a —► 0, where the Reynolds number R has been assumed to be large enough so that the 
coefficients {u m ,Pm} depend only parametrically on the slow streamwise variable x 2 , i.e. x 2 
plays the role of a constant. Substituting (2.2), (1.24) and (1.25) into (1.23) and equating 
like powers of a leads to 


V r - 


‘ VtPo ‘ 
.(tfo-co) 2 . 


<*oPo 

(Uo - Co) 2 


(2.3) 


and 


v T . 


VtPi 

.( Uo-c Q )\ 


a oPi a oPo 0 V r c/ 0 -V r p 0 

(Uo-c 0 y~ zai (Uo-c 0 y + ZCl ( u Q -c 0 y ’ 


where c 0 = So/ao, ai = — iA'/A, and ci = (5i — aiCo)/ctQ. 


(2.4) 


Equations (2.3) and (2.4) must, of course, be solved numerically subject to the bound- 
ary conditions (1.28) and (1.29). However, for the present analysis, it is only necessary to 
know the behavior of the solutions near the critical level. This is most easily determined 
by first expressing (2.3) and (2.4) in orthogonal curvilinear coordinates, say (?;,£), with one 
set of coordinate surfaces corresponding to surfaces of constant base-flow velocity Uo - as 
was done, for example, by Goldstein (1976, pp. 6-10). The functions rj and £ of y and z are 
chosen so that 


U 0 = Uo(x 2 ,r}), 


(2.5) 


rj - y 0 at y = y 0 , r\ -> oo as y -» oo, 


( 2 . 6 ) 
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and 


V r *7 0 -V r C = 0, 


(2.7) 


C = 0 at z = 0, ( = 2n/p at z — 2n/(3. 


( 2 . 8 ) 


where 


yo = < 


0 ; boundary layer 


(2.9) 


— oo ; free-shear layer 

(3 is the (non-dimensional) spanwise wavenumber of the base flow and (2.8) requires (without 
loss of generality) that z = 0 and z = 2n/(3 be the planes of symmetry of Uo. In terms of 
7} and (, the gradient operator in the transverse plane is 


„ Id 

v ' = W”W 


( 2 . 10 ) 


where (Z, m) = (g'Vrj, hV£) are the unit vectors and ( g , h ) = (i V?y| -1 , j V^]" 1 ) are the scale 
factors corresponding to the coordinates respectively. 

It follows from (2.5) that the critical-level position is given by 77 = g c where 


Uo(x 2j g) = c 0 at rj = rj c . (2.11) 

The near-critical-level expansions ofpo nnd Pi can now be found by the method of Frobenius 
(Hall & Horseman 1991; Horseman 1991; and Hall & Smith 1991). To the required level of 
approximation, these expansions are 


Po = aoo + a 02 (7? - 7] c ) 2 + ( 00 ? In |t? - t? c | + b% 3 )(r} - rj c ) 3 

+ (o^ In 1 7? - 77 c | -f a 0 4 + ^X 7 ? “ Vc ) 4 + 0[(r) - t ? c ) 5 In \rj - 7? c |], (2.12) 

and 


Pi = oio + <*11(77 - ?7 C ) + (d[ In 1 77 - 77c | + oi2 + df 2 )(rj - t ? c ) 2 
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+ [(“i^ + d[ L J) In \tj - 7] c \ + 6f 3 ](7? - t ] c ) 3 + 0[{rj - t] c ) 4 In I77 - 7? c |], (2.13) 


where the ± superscript denotes differing values for 77 ^ rj c and the fact that the pressure is 
continuous across the critical layer to O(cre) (see (3.13), (3.17) and (B 3) below) has been 
used. At this point, the coefficients a m 0 and &^ 3 are arbitrary functions of (. Expressions 
for the remaining coefficients in terms of these functions are given in appendix A. 

The boundary-value problem (1.28), (1.29), and (2.4) only possesses solutions for cer- 
tain values of ai since p 0 is a homogeneous solution to (2.3). These values can be found 
without explicitly solving for pi by integrating the difference between p 0 times (2.4) and 
Pi times (2.3) over the transverse domain, applying the divergence theorem to the simply 
connected regions and then making use of (1.28), (1.29), the z periodicity of p and the ex- 
pansions (2.12) and (2.13) to arrive at a solvability condition. For definiteness, we consider 
the simplest case where the critical level forms a single closed or open curve that divides the 
transverse domain into two simply connect regions. In this case, the solvability condition 
becomes 



where the functions and $1 of £ are given by (A 26), the c subscript denotes evaluation 


at T) = pc, 


r2n/p [co n,2-2 

Ip ^L L iv^ 9hd ^ 


Jp 


_ r 2n /P f°° c 0 f \ 

Jo 4 Wo - Co) 3 1' VtPo • v ^° + a ° g °) ghd,ldC • 


and -j- denotes the Cauchy principal value. 


(2.15) 

(2.16) 


9 



For purposes of analyzing the nonlinear flow within the critical layer, it is convenient 


to express the velocity perturbation as 


, . u t v w 

u = z — — + l— -f- m — . 
gh h g 


(2.17) 


The near-critical-level expansions of the shape functions corresponding to u , v and w are 
given in appendix A where it is shown that the discontinuities in (2.12) and (2.13) lead to 
a jump in the streamwise velocity component 


A u = 


30 q 

Cto 


^03 3 * * 6 03 + ^13 ^13) “ ( — — + — l (6 q3 “ &03) 

\O0jj c S>0 «0 / 


+ 


(2.18) 


across the critical layer. Matching this jump with the one induced by the flow in the critical 
layer determines the functions b^ 3 . However, when determining bf 3 , it is more convenient 
to express the jump condition as 


A 


bJonn 


' 02 


u. 


V le 02 +2e u , Uo 


°VV c 


e oi (v - Vc) 


Vc 


(*0 + 3 - %) + *(*>£> - b r 3 ) + (v%- - 21 ) (6& - 

V u o n c 9c a oJ 


= i30 o 

which follows directly from (A 29) and (A 30) 


+ •••, (2.19) 


3. Unsteady flow inside the critical layer 


As already noted, nonlinear effects first come into play locally within the so-called 

critical layer once the deviation of the local Strouhal number from its neutral value becomes 
sufficiently small. The thickness of the critical layer, which is determined by the balance of 
wave-growth and base-flow-convection effects, turns out to be order a on the 77 scale so the 
appropriate scaled coordinate for this region is 

V = 0 ? ~ Vc)/<r- (3.1) 
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The nonlinear terms in (1.9) produce a critical-layer velocity jump at the same order as the 
linear-growth effects when the scale of the frequency deviation cr, which was introduced in 
(1.18), is chosen to be 

<r = c$ ( 3 . 2 ) 

(Goldstein & Choi 1989). Viscous effects will enter into the dominant balance for the critical- 
layer while making only insignificant modifications to the outer flow when the Benney- 
Bergeron parameter 

A = 1 /o z R (3.3) 

(Benney & Bergeron 1969) is order one. In the present analysis, A is assumed to be small 
enough so that viscous effects, which may arise from the x 2 dependence of the base-flow 
solution as well as the viscous-diffusions terms in (1.9), are negligible. 

Since the flow inside the critical layer depends on x and t only through the variables 
(1.21) and (1.22), the appropriate governing equations for this region are obtained by ex- 
pressing (1.8) and (1.9) in terms of a?i , X , fj and £. Upon introducing (2.17), these equations 
become 

a 2 u Xl + aa Q u x + v n + aw^ = 0, (3.4) 

and 

Cu + ghUo n ^ + g 2 h 2 \ 
aCv + h 2 pfj = — a 3 
Cw -f- g 2 p{ — —a 3 


(ap Xl + a 0 p x ) = -a 3 ghAT [ — ) , 


h Kr 
o—N 

( 3 4 

\- 9 -H 2 

h/hyi S ^ 
~-W 

9 

\ h j 

1 9 

9 2 


hw\ 

99C -2 

rrrV 

h>C 

t-w 

h \ 

T/ 

h 2 

h 


(3.5) 

(3.6) 

(3.7) 
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where 



(3.8) 

d v d , d 


~dx + o~dfi + w ~d(i' 

(3.9) 


Af = au - — ■ + a 0 u 

OXi 

Introducing (3.1) into the expressions for u, v, w and p obtained from (A 18), (A29)-(A34), 
(2.2), (2.12) and (2.13) and re-expanding the result shows that the unsteady flow in the 
critical layer should expand like 


u = a^uoi- ui + au 2 + - (3.10) 

v = v 0 + <jvi + a 2 v 2 + •••, ( 3 . 11 ) 

W = <7~ 1 Wq + W 1 + aw 2 + • • • , (3.12) 

V = Vo + crpi + a 2 p 2 + • ■ • , (3.13) 


where, in general, the functions u m , v m , w m and p m of xi , X, 77 and £ have an implicit a 
dependence of the form 

(3.14) 

In this region, the known functions Uq, g and h are given by their Taylor series expansions 
about rj = rjc when expressed in terms of 77 . 

Substituting (3.10)-(3.13) into (3.4)-(3.7) and equating like powers of a leads to the 
following set of equations at leading order 

ao«ox + v 0f j + w 0(: = (3.15) 

£ 0^0 + Uon c ^o 4- &og c h c pox — 0, (3.16) 

Po n = 0, (3.17) 
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(3.18) 


£oWo 4- 

where 

n d , tt - „ v d 

Co 55 C °d, V x + ( a ° U ° - Sl ^dx- (3.19) 

The solutions to (3.15)-(3.18) must reduce to the appropriate linear solutions as xi — > — oo, 
they must be periodic in X, and they must match with the outer solutions discussed in §2, 
It follows from (2.12), (2.13) and (A29)-(A34) that the last condition implies that 

{u 0 ,^o,^o,Po} Re ({i/o-i^/ao^ e 00 , fo-i/fj , a 0 o}Ae iX ) , (3.20) 

as V -> ±oo, where the functions e 00 and f 0 . 1 of ( are given in terms of a 00 by (A 35) and 
(A 44). It is easy to show that the appropriate solutions to (3.15)— (3.18) are 

{u 0 ,u 0 ,u>o,po} = Re ({-t/ 0??c /o-i c £, e 00 Ae lX , ia 0 U OT}c fo-iE , aooAe 1 *}) , (3.21) 

where the function E(x\,X,rj) is determined by 


C 0 E = Ae iX 


(3.22) 


together with the condition that E -* 0 as x x -* -oo and that E be periodic in X. Therefore 

1 /•*! 


e = — [ X1 A({y [x+? «- xl)] az 

^0 J— oo 


(3.23) 


where Y ~ ( a 0 U 0Vc 7j - S^/cq. 

The higher-order critical-layer problems are derived in appendix B. There it is shown 
that the relevant solutions to the order-cr problem can be expressed as 


u U} - + Af) + 


Wi = + w\ -f 




g c h c Ui 


0?j 


c/ f) 


U 0 W o 


gJicUo 


W a 


(3.24) 

(3.25) 
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where the linear components u\- and tu{ are given by (B 18)-(B 20) and the nonlinear 

+ + ^ 
components zi{- and w\ must be determined from 


Oi Q C Q u\- = olqUq^wI^ - ( 73 C - 274 ) Re (iAe lX ) Re(LE^), (3.26) 

= (71 - 272 ) Re (Ae lX ) Re(iE^) - 7 3 Re (iAe 1 *) Re(E n ) - 2'y 2 aoUo r] Re(E) 2 , (3.27) 


where the functions 7 n (C) are given as 


{ 71 . 72 , 73 , 74 } = 


9c 

2a 0 h c U 0rtc 



rv 2 n 2 

a 0 a 00£ 


(3.28) 


and it has been assumed (without loss of generality) that <200 is purely real. The solutions 
to (3.26) and (3.27) turn out to be 


oiou{- = -i( 7 j - 7 3 ) c Re(i^ - G) - ( 7 ^ - 27 4 )Re(iGx + G) - 7 2 C Re(< 7 ) 

+ §( 7 i + 73 ) c Re(^T) + j 2c Re(E)Re(E m ), (3.29) 

= 7 iRe(T> - i F) - 7 3 Re(iF) - 27 2 Re(F)Re(ii^), (3.30) 

where the functions F, G and H of xj, X and fj are determined by 


C 0 {F, G,H}= {Ae ix Re(£ s ), Ae' x ME m ), a 0 U 0vc (F x - 12 F)} (3.31) 

together with the condition that {F, G,H} -> 0 as x\ — oo and that { F , G , Hj be periodic 

in X . Therefore 

f = r r (6 - 

J —co oo ^ 

- A(6)e i[2X+f(6+6 " 2xi)] } d&d&, (3.32) 


G = f X ' t (& ~ fi) 2 ^) 

zGq J—ooJ—oo 


+ -4K 1 )e il2X+r(6+6 " 2 * l)1 } dfid&, (3.33) 
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and 


H = ^ /•*> / & ( Zl _ 6)(6 _ (3.34) 

where the asterisk denotes complex conjugation. 

In appendix B, it is shown that the relevant solution to the order- a 2 problem can be 
expressed as 


V2 m ~ V \ rjrj + a o 


wow? / 


9chcUo Vc + [egjhlUL 


VJ 


4 - u\wq ( U^Wp 


XfjT] 


(3.35) 




9ch c U 0ric T \2glhlUSj__ 

where the linear component is given by (B 36) and the nonlinear component satisfies 
(B 40). For purposes of obtaining the evolution equation for j4(;ei), it is only necessary to 
determine the quantity 


1 f 2 « [ 2 ”/0 . y t 

v m = -J J «ooe-' X 5*„dCdX 


(3.36) 


which, as shown in appendix B, is determined by 


Ul m = J 4Re[(2* 1 + k 4 + 3 fc 6 )4 0) - ( fc 3 + 2* 4 + h)H n - (* a + kt)(E^E^\ 
- (*i + \k 3 - 1* 4 )(£( 1 >4 1 )*) w ] + URe[-i(fc 3 + 2*4 + *5)4°’ 

+ i(*i + k 2 )(E^E^) m ] + 1A*[(2* 1 - * 3 + * 4 + 2* 5 )Gf 



- (*i - ktX&Vnff), - (* 2 + 1*3 + l*4)(£W4 1 >) iW ] - l4 w 

(3.37) 

where 

A 

L s c °dx^ + K^O U 0r?c fj - Si), 

(3.38) 


(• ) <m) = 2^jf'' e “ imX (- ) dX < 

(3.39) 
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the real constants k n are given as 


y*2 7T //? h c 

{ki,k 2 ,k 3 ,k 4 ,h} = / -~{7i 72, 7273, 7i, 7 i 73, 7l) d C, 

JO Qc 


(3.40) 


the function is given as 


<t >2 ~ n J 0 ie lJ: ^ e ( J ^)[^3Re(i r >) - (k 3 + Aj 4 )Re(i.F) 


+ i(2*i + ^ 3 )Re(i^)Re(i^)]dX, 


(3.41) 


and the fact that has been used in arriving at (3.37). The solution to (3.37) 

turns out to be 


vlfjfj ~ (2^i 4- k 3 4 - 3^4 + 4k&)Qi + ( 2ki — k 3 — k 4 + 2k^)Q 2 

+ (2*i - k 3 + Ar 4 + 2 k 5 )Q 3 ~ (k 3 + 2& 4 + h)Q 4 - (k 4 + k 2 )Q 5 
— (k\ — k 2 )Qq — (2Ari + k 2 + |&3 — ^ 4)^7 4" (^2 — 2^3 4- 2 ^ 4)^8 


“ (k 2 + \ k 3 + ^4)^9 - 4> 2 fjfj 1 (3.42) 

where the functions Q n (%i,v) are determined by 

HQi,Q 2,Q3,Q4) = \{AGf\ AGf r , A*G%\ 2/lRc(7/, i )}, (3.43) 

HQ^Qe] = \{2ARe(EWE$\, A’(E^E^ n \} (3.44) 

LW7,<?8,Q 9 } = |{il(J5 (1 )4 ,) *)«. A(E^E^) m , A*(£( 1 >4 1 ) )to}, (3.45) 


together with the condition that the Q n — >• 0 as x 4 —* — oo. Explicit expressions for the Q n 
are given in appendix C. 
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4. Amplitude evolution equation 


The velocity jump induced by the flow in the critical layer will now be computed and 
combined with (2.18) and (2.19) in order to determine the functions 6^ 3 . These results will 
then be used in (2.14) to obtain the governing equation for 4(a?i). 

By using the relation 



2 7TCq 
^qUq^ c 


6(Z-x) 


(4.1) 


where 6 denotes the Dirac delta function, one can show from (3.10), (3.21), (3.24), (B 18)- 
(B 20) and (3.29) that 


/ I OQ . . |T 

& n dfj = - — R e [e^Ae ix )-2n-^ l2i . (*, - f)V(0| 2 df + 0(a), (4.2) 

-oo U U Cq J —oo 

which, when combined with (1.19), (2.18) and (A 36), yields 


^03 ^03 ^ ^ ^03^* (4.3) 

In order to match with the X-independent term on the right-hand side of (4.2), a mean- 
flow component must be included in the solution for the perturbation {u,p}. The ‘steady’ 
Rayleigh problem that governs this component outside the critical layer is given in appendix 
D where it is shown that the corresponding streamwise velocity is of the same order of 
magnitude as the instability wave that produced it and further that the slowly varying 
amplitude of this velocity component is given by 

*= (4.4) 

J — OO 

Again using (4.1), one can show from (B 37) together with the definitions of ot\ and ci 
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that 


£ f r P ^ fc* w.a) 

7T/1 ./o J— oo JO Op V Qp Co / 

f2n/P 

+ 3 7r / aoo #o 

JO 


«1 _ CxE/p^ \ (i) (i) 

ao tO 03 13 , 


d£, 


where 


= -in f 
Jo 

r2n/ 

f R ~-in 

Jo qJJq 


• - l' 2n/P a o9c f g 

9c^Jqt) c 


( ijlc _ 2^^ 
V <7 C 9c J 


a oodC> 


J* s -iTr / 2n//? [s f ^ ^ ] ag) - ^a 02 + P 2 a 0 o 


dC, 


1 77c L \ g e ^U°VcJ ” Sc 
9 = gUoy^/h and V 2 is defined in appendix A . It now follows from (2.19), (3.11), 

(3.36) and (4.3) that 

■2tt/0 


/ 


«00 


2ci <P 
Uorj ( 


~(^03 - Kz) - MKz - Ks - i^aW) 


dC 


2 ^c fr «1 . T Cl 

3o 0 


( / ^ +jh S) + 


Using (3.42) and the results of appendix C, one can show that 


/+00 

-OO 


S *« d9= 4 4 L 


1 n f x i ft 3 


■^(*1 1 £3>&M(6M(6M*(6 + £2 - #i)d£ 2 d£3, 


where 


a. f<&» , ,\ 2 

+ «o«oo d£, 


K = (a?! - 6)[^i(^i - &)(& - 6) ~ * 2(21 - 6) 2 - ^3(^1 - f 2 ) 2 ], 

Vi = 4a 2 0 U^ c (k 3 + 2 k 4 + k 5 ) = j‘""’ h. + <**. 

V 2 = -4c$uS ne (k 3 -ks) = -J o (“o«oo) c df, 

- - " * -*«•.<*. + *.) - - C i (fif 1 ) (% + «w.) «. 


(4.5) 

(4.6) 

(4.7) 
(3.35), 


(4.8) 

(4.9) 

(4.10) 

(4.11) 

(4.12) 

(4.13) 
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Combining (4.8) and (4.9) with the solvability condition (2.14) and using the result 


r2n/p 

I ^o( fl OO®i3^ ^03^1o)^C = 

Jo 


-\2n/(3 


1 V C 


3»o K 


:(aOO a 10C — a 00C«10) 


C=o 


leads to the following amplitude-evolution equation 

rxi p£ i 


where 


K = 


11 = 


rin/p p 

I =I P + I R = 

JO Jy { 

r2n/p poo 


iaoJSt 
(J-I)S o’ 

%ciu 0lc (j-iy 

2n//3 poo 


«oPo 


;ghdr)d(, 


(4.14) 


A' = KA + i/i f ' !' 7f(;ciU 1 > 6)-4(fi)-4(f2)A*(f 1 +6-x 1 )d6d6 (4.15) 

J —oo j — CO 


(4.16) 

(4.17) 

(4.18) 

(4.19) 


lyo (tfo-Co ) 2 

rA7X/p t oo Cc\ f \ 

J = Jp + J R = j ^ J ( _ 3 (VtPo^tPo + alPl) ghdr)d( y 

and the 77 integration in (4.18) and (4.19) is performed along a contour in the complex -77 
plane that lies below the singularity at 77 = r) c . 


Appendix A. Near-critical- level expansions 

In this appendix, the near-critical-level expansions of {u m , p m } are determined by first 
expressing (2.3) and (2.4) in terms of 77 and £. The resulting equations are 


and 


Pom ~ Pon + Vpo ~ 0, 


n q 

Pirn ~ Z — ZTPiv + V Pi ~ ZaiApo + 2 c 17Z — — ^p 0jlJ 


V-Vc 


where 


77 = (77 - 7j c ) 


h 

g(U 0 -c 0 y 


g(Uo — cq)‘ 
h 


J TJ 


C 7 ? - Vc¥ 


= £ HniV ~ »7c) n , 

n — 0 


(Al) 


(A 2 ) 


(A3) 
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V{-) 


1<L 

h d( 


--( • ) 

hdC 


) = £»"(• Kn-VcT, 

n = 0 
oo 

A = a 0 g 2 = A n(r] - r) c ) n , 


n—O 


and 


f? = (g- 7] c ) : 


U, 


O77 


= i>2 n n(V~Vc) n - 


(A 4) 
(A 5) 

(A 6) 


(Vo-toY n=0 

Expressions for V n and A n are easily obtained from the Taylor series expansions of V and 
A about T) — tj c „ The first few coefficients in the near-critical-level expansions of II and Q 


are 


Ho = 2, n 1 = ^ 1 n 2 =(^\ 

9 c \ 9 ) 


U°yr)7) c . ^0 rjr} t 


+ 


Vc 


3 U 0vc * 2 Ug 


(A 7) 


V c 


and 


(A 8) 


Q n — _L. 0,-0 Or, - UoT >Wc _ U °VVc 

0 Un ’ 1 U ’ arp 4f/3 » 

U °Vc ° U 0 t) c 4 U 0ti c 

where g = gUo v /h and the c subscript denotes evaluation at rj = g c . 

Substituting (2.12) and (2.13) into (A 1) and (A 2) and equating like powers of 77 — rj c 
leads to 


®m2 — 2 VoQ, m O) 

a m3 = ~ 3 (T^i U m o — 2J7ia m2 ) 5 

M _ 3 77 ■ (L) 

a m4 “ 4 iil<l m35 

«m4 = ~^[7 ) 2a m o + (Vo - 277 2 )a m2 — -Hi a ( ^ -f 5a^], 

*m4 = 


(A 9) 
(A 10) 
(All) 
(A 12) 
(A 13) 


and 


dn ~ — 2cif?o<zo2? 


(A 14) 
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^12^ — ~^ c 1^0 a 03 •> 


(A 15) 

d i2 = -Oi\A Q a QQ - ciQq (a$ + 3&J 3 ) - \(T1 \d u - S$), (A 16) 

d i3 ~ ^OiiAidQQ -f |ci[2i?2tto2 + ^o( a o4 + 4ao4 -f- 46 q 4 )] 

-U(Po - n 2 )dn - n x {£$ + 24)], (A 17) 

where m = 0,1. 

In view of (2.1), the shape functions corresponding to the normalized velocity compo- 
nents introduced in (2.17) should expand like 


{u,v,w} = {u 0 ,v 0 , w 0 } + v{ui,v u wi} H , 


(A 18) 


as a — * 0. Substituting (A 18) into (2.17) and the result together with (1.24), (1.25), (2.2) 
and (2.10) into (1.26) and (1.27) and equating like powers of a leads to 


and 


where 


and 


io'o'Wo + v 0 n + wo<; = 0, (A 19) 

ia 0 ui + iaq^o + v ir) + c = 0, (A 20) 

{u 0 , u) 0 } = — - — { $Po v , Op 0<: }, (A 21) 

V ~Vc 

{til , W \ } = ^ ^ ^ {$piv* ®Pk} - - c i ^ ) {^o,^o}, (A 22) 

h °° 

*^ (V - ^ c 0 g(U 0 - co) = £ *" (l? ~ rhT ' (A23) 

OO 

6 = - %) -^iWo - co) = £ Gn(n - Vc)n ' (A 24) 

-t OO 

9 = (v-Vc ) 77 - = £ Vniv-Vc)". (A 25) 

U ° ~ c ° iS 
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The first few coefficients in the near-critical-level expansions of 0 and are 


and 


<Z> n = 


= 


&09c 


£>1 = 


U°T)r)i 

W 0r}i 


<*o 9c 

Vo™. Q-n-n . Q n _Ul 


9y c 

9c 


<*o9c V 3C/ o 


Vc 


—2 > 

'Wc __ 9jWc _ Qr] c Uo Wnc , 9yc 

4 Uo„„ % + 9l , 


Vc 


&> = 


aoh c 


01 = 


tfl 




h 


OCqHc 


W 0r}( 


Vc 


hr 


02 


Ut 


OriVVi 


TJ 2 


VV C 


h 


VVc 


&oh c \ SUqjj 4Uq 


y c 


% = 




^ = 






2 C /2 


2 A 


^2 = 


— — 2 > 

hri c Uor)r}c , ^Vc 

» — o 


2 A.CA 


cVQr), 


K 


U 0 


TJ 2 

vvvc 0 nvc 


« 0 J, 


4«? 


*?C 


(A 26) 


(A 27) 


(A 28) 


where h = hUo n /g. It turns out that v m and w m expand like 


vo = e 00 + (4i } ln I 7 / “ »7c| + Cqi)(» 7 ~ 7/c) + (4? ln I 7 ? ~ *7c| + ef 2 )(?7 - 77 c ) 2 

+ ^K 7 ? — ?7c) 3 ln |?7 — ?7 C |], (A 29) 

v x = e$ ln | ?? - Vc \ + ef 0 + (eg> ln |t? - »y c | + e*)^ - ?? c ) + 0[(»y - ?? c ) 2 In \g - Vc \], (A 30) 
and 

= fo~i(v - 9c)~ l + foo + foi(v - 9c) + 0[(tj - 77c) 2 In |t? - ifc)], (A 31) 

= fi-2(9 ~ 9c)~ 2 + /i-i(t? “ 9c)~ l + ho + 0[(g- g c ) In |7? - 77 c |], (A 32) 


as 77 — ► 77 c, where the coefficients e mn and f mn are at most functions of xi and (, It therefore 
follows from (A 19) and (A 20) that 


uq — i« 0 1 [/o-ic( 7 7 ~ 9c) 1 + e oi In {9 — 9 C \ -f + ej^ 4- foo ^ 

+ ( 2 e 02 ) ln I 7 ? ~ 9c\ + £02 + 2^02 + foic)(9 - 9c)] + 0[(9 - 9 C ) 2 ln (77 - g c \], (A 33) 
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1 [/i-2c(»7 “ Vc)~ 2 + (4o + fi-H ~ a i a o Vo-ic)(»7 “ Vc ) _1 

+ ( e n “ ai«o le m } ) ln I 7 ? ~ Vc\ + eff + e± + f 10< - + e± + / 0 o c )] 

+ 0[(v~ 7/ c )ln|7?- ?? c |], (A 34) 

as 77 -+ Tfc. 

Substituting (2.12), (2.13) and (A29)-(A32) into (A 21) and (A 22) and equating like 


powers of 77 — tj c leads to 

eoo = i2<I>o a Q2 , (A 35) 

e oi^ = (A 36) 

e oi = i[2<Mo2 + ^o(«S)3 + 3 6±)], (A37) 

e 02^ = i(3 ^i«03^ + 4 $o fl o4^)? (A 38) 

e 02 = ^[2^2^02 + <Pi(a\g + 3 6q 3 ) + + ^ a o4 + 46q 4 )], (A 39) 

4o = i2 ^0^12 + ^1 ^4? » (A 40) 


efo — 4- ^ 0(^12 4* 2ai2 4- 2^)] — aiao le oo + Ci(^ieoo 4- ^o e oi)> (A41) 

4? = i[2^4a + 30b(ag } + ^3)] ~ *4? + c i( *i4? + ®b4a), (A 42) 

e n = i[^2^n 4* ^i(^W + 2«i2 4- 2^2) 4- $o( a i3^ 4- 4.3* + 4-36j^)] 


— oiOq le oi 4- Ci( ^2^00 4- ^i^oi 4- $o e o2)> (A 43) 

fo-i = i 0oaoo^ 5 (A 44) 

/oo = i©i a 00 £, (A 45) 

/01 = i(0 2 «ooc 4- ©qo,q 2 ()’> (A 46) 
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and 


/ 1-2 = Ci^d/o-i, (A 47) 

fi-i — i<9o a io^ “ ^ 1^0 i/o-i + c i(^i/o-i + ^o/oo) 5 (A 48) 

/io = i(^i a ioc ^o^nc) ~ a i a o 1 /oo + Ci( $ 2 / 0-1 + ^ 1 / 00 + S'o/oi)* (A 49) 

Appendix B. Higher-order critical- layer problems 


In this appendix, the higher-order critical-layer problems obtained by substituting 
(3.10)-(3.13) into (3.4)— (3.7) and equating like powers of <j are given. The order-<r problem 
reads 


a 0 u lX + ^0^! + Vifj -f Wi^ = 0, 

£oU! + CxUq + u 0r)c + h^nv oj + 9ch c {otoPix + Po Xl ) + 2a 0 (gh) Vc fjp 0x = 

Pirj = 0 , 

C0W1 + £iw 0 -f- y P i < + 2 ^^°C = ~^i» 


(Bl) 
-0i, (B2) 
(B3) 
(B4) 


where / = ghU 0 n , 


r _ , TT - ^ , 1 rr _2 ^ 

1 - 0 + + 2“ ol/o ^ ,? ax’ 


-e- 

H-i 

III 

P 

o 

f u o > 
\9ch c ) 

K + 

/ U 0 U 0 > 

V gchc j 

>.♦ 

{ u 0 u>o\ 
V 9ch c J 

o 

a 

III 

qT 

'u 0 w 0 \ 
- 9ch c J 


'v Q Wq' 
k 9ch c j 

)_ + 
V 

iN' 

h c y 9c ) 


It follows directly from (B 3) and matching with the outer linear solution that 


(B5) 

(B6) 

(B7) 


Pi = Re (a ]0 j4e lX ) 


(B8) 
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It turns out that, for purposes of computing the velocity jump A u across the critical layer, 
it is only necessary to know u lf} , v li}fj and w l . Therefore (B 1 ), (B 2 ) and (B 4 ) are rewritten 


as 


- u\ n ) x + v lm - v\ m + (uq - u)})_ c = 0, (B 9) 

C Q {u ln - u}_) = 0b„ c (®i - tDl) c - (B 10) 

£o(wi - w\) = -6 U (B 11) 

where u}-, v\~~ and u)} satisfy the linear equations 

a ou\ nX + u 0xin + v\.. + w\- c = 0, (B 12) 

+ (f;i7 p 0<) + “ 0 9ch c (^ - 2^) Poxx = 0, (B 13) 

£ o »l + /:iWo + | £ Pi ( + 2^poc = 0, (B 14) 

and have the following large - 77 behavior 

~ Re ({ieoi’/oo^, 4n/v, fm}Ae ix ) (B 15) 


which ensures that the solutions to (B9)-(B11) match with the outer linear solution as 
fj -*■ ±oo. By using (3.21) and (3.22) together with the relation 


A(* ) = C 0 [M 1 (- )]- 


UPVc _ ^ 0r )V c \ 

Co 2 £ 1 0 „J n So 


4-Co( • ), 


(B 16) 


where 

M(-) = 

it can be shown that 


Uorjc Uor, nc \ _ Si 


c Q 2 U Qvl 


So 




(B 17) 


a 0 u \fi — ( u 0 X 1V + jjfj + w \f}0x^ 


(B 18) 
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v i fjfj = ReiiaoUo^etfE), 


i — — ^io:o^07j c /o-i M.\ + 2^^-aoo^r) -f j^-a E 


It follows from (B 6), (B 7) and (3.15)-(3.18) that 


, tto _ ( U 0 W 0 \ 

r’l t] Tj POX^Orir) "b ( L ) ^0 

u 0r)c \g c hcJfj<: 


/ U 0 U 0 f j\ 

\ f c J 


VJ 


1 «0 1 I w 0\ 

1 - w^ c po < uo * ~ cw Pox “°' + (r L _ £o 


U 0 w 0 

fc 


(B 19) 


(B 20) 


(B 21) 


(B 22) 


Combining these expressions with (BIO), (B 11), (3.21), (A35) and (A 44) then leads to 
(3o26) and (3.27). 

The order-(7 2 critical-layer problem reads 


<*oU2X + Ui Xl + v 2f j + w 2 c = 0, 


(B 23) 


£oU 2 + Citti + C 2 U 0 + U 0 TI C («2 + -jr r l v l + + Sch c (aop 2X + Pi*,) 

+ 2(ffft)„ c ^(«oPix +P0x,) + 


{ ^Q^^ aoPox = (B24) 


t C °* o + ^ = ^o 


h r 

dlhr 0 ' 
9 d c ^ , to 2 )? 


(B 25) 


£o»2 + £ 1 % + £ 2 «>o + | £ P2< + 2^#1< + ^nr V 2 Po C = -02, (B 26) 


where 


Co = 


to*0 


»7»7c 


We 


2s./»c 


(gh)Ti C:= ~ , l 


0 l 


0 


*/ 2 £°+ gchc Cr lCi + 2 Uo, ”'^ 2 g Xl + QVoUonrmJ? gx , (B27) 


^2 = 2«0 f 5°5i') + /'Sl®0 + »0®l\ | f Slgp + tjggA 

\9ch C J X V g C h C J fj \ gjlc J £ 


+ 


9chc 


Xl 


(g h )v c - - 1 r(Acl ~ - 


(B 28) 
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fl 2 __ ao / «qwi + ^!^o \ ^ / uqu?i + ^iwq \ 2 

\ fl^c / jf \ Pc^c / 77 he \ 9c ) £ 
fuowo\ _±(l) y niDn _J__.fiVc\ 

9c\h) Vc 9chc\9c) t 


Wo- 


(B 29) 


\ 9ch c V a;! i/c \"V ? 7 C i/c"c \ ye / ( 

Fortunately, only the solution for v 2 ^ is needed in determining the governing equation for 
4(ari). Therefore the above equations are combined to give 


£0(^27777 — ^27777) — [ttofox + i’lxx + #2< + — wl)]jj 


, T 1 t' 


, (gh) - a fr 


d 


P v I /V c f* ^ ^ rr 

£1773? H r — — a o~T^^On r -^77 

or} 2 g c h c dg f "*dX 


{di ~ 4 )+ °{i£) (bso) 


where 


and vl^ is determined by the linear equation 

^°^2t777 ^17777 ^ “ /, ^ T £ ^ r 0r7 c ( Q; 0^ix + ^0a:i ) ~ (AW* + £2^0)77 

J C 


y _ ® (9c d ^ . o . d 2 




9c O \ 2 

fc c 0cJ + “° 5cAc 


ax 2 ’ 


(B31) 


9chc 


f m c 


^° v °) - r/ <j aa'-j^-Uo^vox + 


(<7 2 ) 


m c 


P0( 


, 2 (0 2/ * 2 W„ 

+ o n ; Poxx 


9ch c 


= 0, 


+ £>., _. - ■ ■ g chc 

— 4 a o(gh)Ti c POx 1 X “ 2 + Q; o(5 f ^)r7 c Pl J YX 

together with the boundary condition 

^2^77 ~ * ^ e [(2eQ 2 ) In \afj\ + Se^ + 2e^ 2 4 /fj)Ae lX ] as fj — > ±oo 


(B 32) 


(B 33) 


which ensures that the solution to (B30) matches with the outer linear solution. By ma- 
nipulating (3.15)-(3.18) and (B 1)-(B4), one can show that 


Uorj'.i&ov \x 4 — £0^1^ — g ( a o~r^UoT] c v ox H ^iPo + 2— ^£>oPo ) 

\ J c 9c 9c J 


- 2o:og c h c po Xl x - yV 0 pi 


(B 34) 
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and 


(Ac™! + C 2 ^w 0 )- = - 


( 9 h ) 


Vc 


9cK 


2 V 


(9 


Vc 


Poc + ~P H 


~ V 


(gh) 


Wc 


L 9ch c 


TPo< (B35) 

£ Me 


where the V n are defined in appendix A. Combining these results with (3.21), (B 8), (B 16) 
and (B 19) then leads to 


V 2r)r) 


V 9 C he J 


.-,t 


U 0 


v[ n + ~^v\ri - Re (vi’oV 0 a 02 Aj x ) 


- Re ( 2 ^ [too - 2&) gla 0 0 A + («1 E ) . ( B 36) 

1 9c L \ 9c 9c J Cfx i g c Z J J 

For purposes of computing the induced velocity jump, it is convenient to express (B 36) as 




Ut 


° T ^-v\~ — Re 


9 V2m U 0nc ''*« 


J. L ) i o p ± _ 

c 02 ' ■ £C n9 


U, 




02 tt C 01 
U °V C . 


Ae 


iX 


= — Re < 3 


A 

9c L 


„(£) A (l) , : 2 / *"/ c 9 

Utj — «Q3 +177^0 1 ^ 


W ( Ac o Ac \ „2 

13 “03 ' 1 ^ 




, ^c^OO^- | 

9c J dx x \ 


E 


- Re { 2— 
l 9c 

+ ReiA<i (1) 


- &) a\S + - » 2 a 00 

2 c o 9c * 9c 


Lfij 


Sc 


03 


Uqtj c 


[\ c 0 2 u 0rlc 


f * 


L Er 


(B 37) 


where 


‘‘1 5fc( a -0' 

and (3.22) and (B 19) were used in arriving at (B 37). 

It follows from (B28), (B 29), (3.15)-(3.18) and (B9)-(B11) that 

{<Xoi>2X + - u, — ^ hi” j 


(B38) 


U 0 


(pox*i w ), 


+ 


_1_ f 2 wow i ^o^g 
A \ 9c 9cfc 


,t,2' 




+ («0^2X + ^lc) 


£o i 

L 

u 0 u\ , 

r + 

f4' 

)' 

+ 


( “o^o \ 

► 

1 

l 

_ fc ' 

wl, 

1 - 
VJ 


fc 

l 2/1 

c. 


(B 39) 


rm 
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where V >2 and 0\ are given by the right-hand sides of (B 28) and (B 29), respectively, but 
with {u lt v u wi} replaced by {u}, v\, u){}. Introducing the above relation into (B30) leads 


to 


nn ~ ~y7^~(Pox u lfjfi) y ~ 


ct. 




Uo Vc ^ lm, x U 0r)c 


vv 


<*0 


(p ox w 1 TO ) ( 


hi 


Ot] 


+ 


1 (* w o™\ . u 0fj wl' 


K l 2 9c + g c f, 




+ [ a 0^h + VOan + + A((wi — U)|)]_ 


r 1 L , 

c\~o T 


Co-h - oco^U 0 „~ 


dv 2 ' 9ch c C °df) U ° f r U °^dX 


(Sl -* 1) + e fe?) (B40) 


where v\ m is given by (3.35). Substituting (3.21), (3.29) and (3.30) into (B 40), multiplying 
the result by a ao e~' x / it, and integrating from f = 0 to “in /f) then using the relations 

[2a 0 {/ 0 „ c Re(£)Re(fy)] Xs = Re (rle iJr ) [Re(Fy) + Re(£)Re(i^)] w 

+ Re (iAe 1 *) [Ke(E)Re(E„)] m - £ 0 [Re(£)Re(f Jf )] w , (B 41) 


[2a 0 Efo, I< .Re(.E)Re(if’)].yjj = Re (Ae 1 *) [Re(if)] w 

+ Re (iAe 1 *) [Re(£)Re( £„)]„, - £ 0 [Re(£)Re(iF)] w , (B 42) 

[ a o b'o, c Re( 1?) 2 Re ( i E. Fl )]xxij = Re (Ae lA ) [Re(i£)Re(i£^)] w 

+ Re (iAe ! *) [Re(B)Re(i£,-)] w - £o[Re(£)Re(i£)Re(iR; fl )],„ (B 43) 
and integrating from X - 0 to 2 n leads to (3.37). 


Appendix C. Expressions for the Q n 

The solutions to (3.43)~(3.45) are 

M f x i ft 3 „ 

<?i = -*"5- / / / (f2-6) 3 Ci(x I ,o|&>£2,6)d{idkd£ 3 , (Cl) 

** J — CO J — oo J — oo 
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M t Xl f^ 3 r, 

<32 = 1-5-/ / / (f 2 -6) 3 C 2 (2r 1 ,i;|af2,6)dfidf2d?3, (C2) 

6 J—oo J—oo J—oo 

. M /* Xl ft 2 

Q 3 = 1 ~o' / / (6 - 6) 2 (2& - 6 - 6)CMzi,^|f3>6>6)d6d6d6, (0 3) 

^ J—oo J—oo J—oo 

Qi = iM f I f (6 - 6)(6 - 6) 2 [C'i(*i»vk3»6»6) 

J—oo J —00 J—oo 

- C 2 (*i,i?l&,f2,fi)]dfidf 2 d{ 3 , (C4) 

/ ^i /*^3 r^3 

/ / (6-fl)[(&-f2) 2 

-00 7 — 00 J — oo 

+ (?3 - fi) 2 ]C'i(£i,i/|f3>f2>{i)d{id{ 2 d{3, (C5) 

Qe = iM f X ' f /' £3 (f 3 -{i) 2 (2f3-&-f 1 )C 3 (a:i ) 5|f3,f 2 ,f 1 )df 1 df 2 df3, (C 6) 

J — oo 3 -oo J—oo 

Qr=-iM P f f (6-6X6 - tifc^xun I&.6. fi)dfidf 2 df 3 , (C7) 
Qs = >M r l / (f 3 -6)(6-fi) 2 C’i(x 1 ^|&,6,6)dfid6df3, (C 8) 

J—OO J — OO J—OO 

and 

<?9 = i M [ 1 I P (fo - 6)(26 - 6 - fi) 2 C 3 (*i, 5 1 f 3 , 6, fi)d{idf 2 d&, (C 9) 

J—oo J—oo J—oo 

where 

Ci s ^(fe)2l(6)A*(f 1 )e i? < &+?2 - ?1 - I1 >, (CIO) 

C 2 = ^(6)4*(6)^(6)e iFK3 " <2+fl ’ xl) , (Cll) 

C 3 = A*(6)4(f 2 ) 2 l(fj)e ip( " fo+{2+fl - I,) , (C 12) 

and ilf = o,lV^j2c%. 

By using (4.1), one can show that 

/ +°° r x i rt 3 , 

Qidr] = —IN / / - £3) -P(^i | £3 ? £2)d£2d£3, (C 13) 

-00 J —00 J —00 
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f +° o f+°° rx i riz 

/ Q 4iy= Q a dr) = iN / (si - f 3 )(& - 6)£(*i 1 6,6)d6d6. (C 14) 

J—oo J —oo J—oo J—oo 

/ +°° m /*6 

Q 5 dn = -iN / (x 1 -&)[(x 1 -6) 2 + (6-6) 2 ]£ , (^iU3.6)d6d&, (C15) 

-oo J—oo J — oo 

/ +°° [X 1 

(Mi? = -UV / / (xj - f 3 ) 2 (zi - 6)£>(*1 1 6, 6)df 2 d&, (C 16) 

-oo J —oo J —oo 

and 

/ +oo r+oo r+oo r+oo 

Q 2 df]- / Q 3 dfj= / / Q 9 d?? = 0, (C 17) 

-oo J—OO J—oo J—oo 

where 

-D = A(^3)A(^ 2 )^*(£3 + £2 - £ 1 ), (C 18) 

and N = n agt^/cg. 

Appendix D. Mean-flow distortion 

In this appendix, the solution for the mean-flow distortion generated by the critical- 
layer nonlinearity is analyzed. When the mean-flow distortion terms are made explicit in 
(1.19) and (1.20), these equations become 

u = Re ^AuQ iX j -f Re ^iB—^ -f IcrB'j^ -f maB'— ^ + . . . , (D 1) 

p = Re ^Ape 1 ^ -f Re(cr 2 B"p) -f . . . , (D 2) 

where B(x i) is a slowly varying amplitude function and the functions u, v , w and p of ®i, 
y and z expand like 

{u,v,w,p} = {u 0 ,vo,w 0 ,po}(y,z) + •••, (D3) 

as a 0. Substituting (D1)-(D3) into (1.15)— (1.17) shows that po satisfies the ‘steady’ 
Rayleigh equation 

v -Rf)=°- (° 4 ) 
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while the velocity fluctuations are determined in terms of po by 




h „ 9 - \ 

Pot), ~J L P°Cj 


(D 5) 


Near the critical level, po expands like 


Po = r 00 + r^(rj - T) c ) + . • • 


(D 6) 


where (3.21), (B 8) and (B 25) have been used to conclude that the mean pressure fluctuation 
is continuous across the critical layer to 0(a 2 e). It follows from (D5) that the discontinuity 
in (D 6) leads to a jump in the streamwise velocity component 


Au 0 


hcUor) c 

9ccl 


( r oi — r 0l) 


across the critical layer. Matching this jump with (4.2) yields 


(D7) 


01 


01 


= —2 71 


9cl2C 

h c CQ 


and the amplitude equation (4.4). 


(D8) 
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